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ABSTRACT 


A  model  for  the  interface  of  two  interacting  one  component 
plasmas  of  different  background  densities  separated  by  an  impermeable 
surface, is  solved  exactly  in  two  dimensions  at  one  special  temperature. 
This  model  emulates  the  behaviour  of  an  ideally  polarizable  interface. 
We  obtain  density  prof iles ,d iff erential  capacity  and  surface  tension  as 
a  function  of  the  applied  potential. Pertinent  sum  rules  are  given. 


— * 


I  INTRODUCTION 


A  model  for  a  charged  interface  has  been  proposed  recently  (1,2) 
in  which  two  media  described  by  classical,  one  component  plasmas  (OCP)  of 
different  background  densities  are  separated  by  a  charged  gap.  Although  this 
model  is  far  from  a  realistic  description  of  the  metal-electrolyte  interface, 
it  has  the  interesting  feature  that  an  exact  solution  can  be  found  in  two 
dimensions  for  the  special  reduced  temperature  of  2, using  the  method  of 
Jancovici  and  coworkers  [3] .On  the  other  hand,  as  shown 
in  [1]  ,it  appears  that  this  model  does  not  allow  any  change  of 
the  potential  drop  across  the  interface  with  the  surface  charge,  and  in  this 
sense  emulates  the  behaviour  of  a  non-polar izable  interface. 

This  is  a  result  of  the  fact  that  charged  particles  can  move  freely  across 
the  interface, which, because  of  perfect  screening, is  probably  connected  to  the 
non-existence  of  long  range  correlations  along  the  interface  [1,2], 

In  the  present  note  we  consider  another  physically  interesting  situation 
corresponding  to  the  polarizable  interface,  such  as  the 
classical  example  of  the  mercury-aqueous  electrolyte  solution  interface. 

For  that  purpose, we  assume  that  between  the  two  media  there  is  an  impermeable 
surface  which  prevents  the  ions  from  crossing  the  interface  while  allowing 
electrical  interactions  between  all  particles. Clearly,  this  model  is  consistent 
with  the  concept  of  the  ideally  polarizable  electrode  [4]  .A  full  description 
of  the  technical  details  of  our  solution  will  be  given  in  a  forthcoming 
paper. Here  we  give  only  a  summary  of  our  main  results. 


II-  PARTITION  FUNCTION:  METHOD  OF  SOLUTION 
Our  system  consists  of  two  concentric  circles  of  radii  R^  and  Rj 
The  inner  region  (0<r<R^)  has  background  charge  density  -eoc.^  , where 
e  is  the  elementary  charge,and  contains  N^  partic’  les  of  charge  e 
We  want  to  allow  for  an  excess  charge  density  -eo*  along  the  wall 
then  we  must  have 

K,-  (i) 

The  outer  region  (R^ <r<R^ )  has  background  density  -eof^  rexcess 
charge  density  eo^  and  contains  N^  partic.  les  of  charge  e. Hence 

^  (2i 

The  plasma  coupling  parameter  is 

r'=/i«t=2 

where  (i*l/kT  is  the  thermal  Boltzmann  factor  and  T  is  the  temperature 
The  potential  energy  of  the  system  is  ^  .  ^ 

\J  -  \]q  +  J-ZZ  f*  I?-*}  I *■  +  ZZ  ^  +  ZZ 
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where  V0  is  a  constant  background  term  and  N^N^+N^  .Note  that  the  charge 
density  eo'  does  not  appear  explicitly. 

Following  the  method  introduced  by  Jancovici  [3] ,the  canonical  partition 
function  can  be  written 
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where  is 


the  permutation  operator  between  the  N  particles. 
Since  the  ions  are  not  allowed  to  cross  the  inner  circle, it  will  be 
useful  to  introduce  a  factor  z  which  ,  as  we  will  see  later 
fixes  the  ratio  of  the  activities  on  both  sides  of  the  interface, 
can  then  be  written  as  a  contour  integral 
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with  hz  s  n  + 

and  z  J  n*  Ju. 
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is  the  incomplete  gamma  function.lt  is  convenient  to  write  (5)  in  the  form 
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"i  and  H  are  the  partition  functions  for  the  uncoupled  inner  and  outer 
regions, and  are  introduced  to  insure  proper  convergence  of  the  integrals 


below 
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Define  K,  »< 
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In  the  limit  R_  ,R.  ->  oo,the  interface  becomes  a  straight  line,  and  we  get 

y<^~  R,  ^  I  /  V*.  r±jJS!^ljb  +  [  -fit  f  jfc  ] 
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In  this  limit  ^  can  be  evaluated  asymptotically  by  the  method  of  steepest 
descents  , choosing  a  contour  z*z#exp(i9).  We  find 
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where  z#  is  the  position  of  the  saddle  point  in  the  complex  z-plane. 


z *  is  real  and  is  computed  from 

,0o 
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For  m*l  we  have 

*.(-*0*  *0c*y 

so  that  Z#(«-«o\s  1 

In  the  general  case,  (10)  must  be  solved  numerically. 


(10) 


III-DENSITY  PROFILE 


The  density  pit)  can  be  obtained  by  the  same  method. After  some 
manipulations  we  get 
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for  r<R1  (11) 
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for  r>R  (12) 
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In  the  limit  R  ,  R  ->oo  and  taking  the  origin  of  coordinates  at  the  wall 

X  ^ 

(xs  r-R  )  ,we  get 
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The  integrals  in  (13)  can  be  computed  numer ically .The  results  for  one  typical 
case  are  shown  in  Figure  1.  For  m=l  the  contact  densities  have  the 
simple  expressions 


(14) 


Unlike  the  case  of  the  permeable  [1,2]  inter  face , the  density  has  now 
a  jump  at  the  interface. Notice  also  that  for  zero  charge  and  m*l,the 


density  remains  constant  everywhere. 


IV-POTENTIAL  DEOP  AND  DIFFERENTIAL  CAPACITY 


The  total  potential  drop  across  the  interface  can  ae 

calculated  by  integration  of  the  density  profile,  however,  this  is  not 
necessary,  since  b±ip  is  directly  connected  to  which  is  the  ratio 
of  the  activities. Indeed  we  have 


-  m  fi 


-a 


where  and  are  the  chemical  potentials  on  each  side. In  fact 
it  appears  that  is  the  true  external  variable  controlling  the 

state  of  the  system. This  is  a  very  nice  feature  of  this  model, 
because  this  is  actually  what  happens  in  the  electrochemical  measurements 
where  the  potential  is  externally  fixed. For  a  given  bxifi /  zo  is  obtained 
from  (15),  and  the  charge  density  from  (10) .Alternatively , using  the 
electroneutrality  condition, we  have  a  simpler  way  of  calculating  ecr 
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Eq.(16)  can  be  shown  to  be  completely  equivalent  to  (10). 
The  differential  capacity  o  is  now  obtains 
form  from  (16)  ^ 


is  now  obtained  in  closed 
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A  typical  curve  showing  the  dependence  of  C  with  potential  is  shown 
in  Figure  2. 

For  m-1,  C  sl{% 


V-THERMODYNAMICS  AND  SUM  RULES 


From  (6)  and  (8)  the  surface  excess  Helmholtz  free  energy 


can  be  calculated.  We  get 
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where  f*  is  the  surface  excess  free  energy  of  the  uncoupled  system 
Differentiating  with  respect  to  the  surface  charge  density, we  ind 


the  first  sum  rule 

+  ^  ^  U9) 

It  may  be  shown  that  this  equation  is  equivalent  to  the  statement 
that  our  model  satisfies  Lippmann's  eguation 

(2z\  -«£. 

yr  7,1  (20) 

where  ^  is  the  surface  tension  defined  by 


-c^) 

^  L(,T,V/ 


where  Q»edA  is  the  surface  charge. 

A  similar  result  holds  for  the  case  of  a  OCP  near  a  charged  plane  [5  ,SJ. 
Figure  3  shews  the  dependence  of  the  surface  tension  with  potential 
(electrocapillarity  curve)  . 

Another  sum  rule  can  be  derived  by  thermodynamic  arguments: 

*?;  2 -lt  [  f,«s\  -  f&S]  + 


where  9.  and  P.  are  the  bulk  thermal  pressures  on  each  side  of  the 

1  l 


inter  face. Notice  that  this  contact  theorem  differs  from  that  of  the  permeable 
interface  case  [7]. Equation  (22)  is  indeed  satisfied  by  our  microscopic 
expressions. 

Additional  sum  rules  and  a  more  comprehensive  discussion 
of  this  model  will  be  given  in  future  work. 
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